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ABSTRACT 

The evolution of compact groups of galaxies may represent one of the few places in the nearby 
universe in which massive galaxies are being forged through a complex set of processes involving tidal 
interaction, ram-pressure stripping, and perhaps finally "dry-mergers" of galaxies stripped of their 
cool gas. Using collisionless iV-body simulations, we propose a possible scenario for the formation of 
one of the best studied compact groups: Stephan's Quintet. We define a serial approach which allows 
us to consider the history of the group as sequence of galaxy-galaxy interactions seen as relatively 
separate events in time, but chained together in such a way as to provide a plausible scenario that 
ends in the current configuration of the galaxies. By covering a large set of parameters, we claim 
that it is very unlikely that both major tidal tails of the group have been created by the interaction 
between the main galaxy and a single intruder. We propose instead a scenario based on two satellites 
orbiting the main disk, plus the recent involvement of an additional interloper, coming from the 
background at high speed. This purely iV-body study of the quintet will provide a parameter-space 
exploration of the basic dynamics of the group that can be used as a basis for a more sophisticated 
iV-body/hydrodynamic study of the group that is necessary to explain the giant shock structure and 
other purely gaseous phenomena observed in both the cold, warm and hot gas in the group. 
Subject headings: galaxies: evolution — galaxies: interactions — galaxies: individual (NGC 7318a, 
NGC 7318b, NGC 7319, NGC 7320c) — methods: numerical 



1. INTRODUCTION 



More than 130 years after its discovery (JStephan 18771 ). 
the Stephan's Quintet (aka HCG 92 or Arp 319; here- 
after SQ, see Figure [lj remains enigmatic, despite be- 
ing one of the most well studied compact groups in 
the local universe. One of the particularities of SQ 
is the wide range of redshifts of its supposed mem- 
bers. The southern galaxy (NGC 7320) has a velocity of 
800 km s~ \ i.e. about 5000 km s -1 les s than the rest of 
the group (Burbid ge fc Burbidgdll96ll) . Although early 
hypotheses argued that this ga laxy was a case of exis- 
tence of non-Doppler redshift (|Arpl Il973al ). it is more 
likely a foreground galaxy, as there are no evidences of 
interaction with other members of the quintet. This 
has been confi rmed by redsh ift -independent distances 
measurements J[Shostaklll974ai rbl: lAllen fe Sullivadll98l 
iMoles et al.l I1997D . and by the clear resolution of indi- 
vidual stars in NGC 732 by the r ecently refurbished 
Hubble Space Telescope (|Nolll[200l ER04). However, 
the higher-redshift system can rightfully be classified as 
a quintet because of the discovery of an additional group 
member, NGC 7320c, which clearly shows optical links 
with NGC 7319 and likely was responsible for at least 
one major t idal filament in the group as suggested by 
lArjl (jl973bf ). Finally, one of the members, NGC 7318b 
appears to have recently fallen into the group from the 
background and is currently interacting with the intra- 
group medium (IGM). 
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The puzzling aspects of Stephan's Quintet do not 
end with its kinematics. It has long been known 
that its cool gas distribution is very peculiar. Almost 
all of the H I emission is seen outside of the opti- 



cal body of the group members (lAllen fc Sullivan] 



1980; 



1999 



IShostak. Allen fc Sullivanl [TMl IWillia'rns' et al 

2002), and even its cold molecular gas is mainly displaced 
from the galaxies and lies in apparently isolated clumps 
in the IGM, with only a few concentration within the 
optical structure of the largest galaxy NGC 7319. 

Even more peculiar was the discovery of a radio con- 
tinuum filament lying be tween the NGC 7 3 18b/a 
pair and NGC 7 3 19 (lAllen fc Hartsuikerl Il972l : 
Ivan der Hulst fc Rotsl 119811 : Ku et all 120031 ). per- 
haps suggesting a shock compression caused by the 
interaction of the infalling galaxy NGC 7318b with 
the group. The filament was su bsequently found to 
be a strong emitter of X-ravs dPietsch. et all 119971 : 
iTrinchieri etaI1l2003L 120051: IP 'Sullivan et al.ll2009h . and 
optical spectroscopy provided more evidences that the 
struc t ure wa s a giant intergalactic shock wave (|Xu et al.l 
l200l 120031) . Surprisingly, the X -ray luminosity in 
the shock structure is exceeded by emission from 
lines of warm molecu l ar hydrogen seen in the Mid-IR 
(|Appleton et al.l 12001 iCluver et all [2010h . All of these 
peculiarities seem to relate to a complex and rich history 
of collisions and interactions in the group which has led 
to a highly disturbed state for both th e gas and stars 
in this system (e.g . IMoles et al.l Il997t iSulentic et al.l 
120011 : iGuillard et all 120091 ) . Attempting to unravel the 
complex past history of the group via simulations is the 
main aim of the present paper. 

Our understanding of the dynamical prop- 
erties of compa ct gr oups in general (e.g. 
Barnei Il985t iDiaferio. Geller fc Ramellal 

19941 : lAthanassoula. Makino. fc Bosmal ~ 119971 : 
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800 km/s 

Fig. 1. — Top: Stephan's Quintet members and the foreground galaxy, observed in the B-band. North is up, East is left. Image courtesy 
of Vicent Martinez, Fernando Ballesteros, Vicent P eris and Rodney Smit h. Bottom: schematic view of the main morphological features of 
the group. Velocities of the galaxies are taken from lSulentic et" al. ( 2001). 



lAceves fc" Vclazqucz 2002]) relies on statisti- 
cal explorat ions covering broad ranges of pa- 
rameters (IMcConnachie, Ellison. & Pattonl 120081: 
iDiaz-Gimenez fc Mamonl I2008D . Although this numer- 
ical and theoretical approach provides an important 
global description for the interaction process, it misses 
a comparison with compact groups actually observed. 
Therefore, while the physics involved in the evolution of 
such objects is well-known, the question of scenarios that 
might bring a new light on observations of individual 
groups remains poorly addressed. 

In this contribution, we investigate the possible chains 
of events that may lead to the observed configuration of 



SQ. Because of the intrinsic challenge that represents the 
simulation of the group, we concentrate first on the col- 
lisionless dynamics, leaving the additional level of com- 
plexity of the hydrodynamics to a forthcoming paper. 

We first draw a quick picture of the group, listing the 
features we use as clues to understand the formation his- 
tory. In Section 3, we propose an original approach to 
explain the web of features of SQ. We distinguish two 
main regions in the quintet and describe the formation 
scenarios used in Section 4. We close with a discussion 
and a summary of our main results. 

2. MORPHOLOGICAL FEATURES 
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Early observations showed tha t H I is widely spread 
outside the optical regions of SQ (|Allen fc Sullivanlll98fl ; 
iShostak. Allen fc Sullivanlll984f) . One classical explana- 
tion is that the gas has been stripped-off the main galax- 
ies of the group during tidal interactions and close pas- 
sages. More recent measures (jWilliams et al.ll200l con- 
firmed that the radial velocities of H I regions are in good 
agreement with this idea, as they match the velocity of 
their former host galaxies. 

Optical images show two large tidal arms, both ex- 
tending from NGC 7319 to the southern region of the 
group (in blue on Figure [T]). The outer one, very dif- 
fuse, passes behind the foreground galaxy (NGC 7320). 
Another arm, more concentrated, is almost parallel to 
the first one and seems to point to the small NGC 7320c 
galaxy. The differences in length and luminosity of these 
two structures have been explained by two separated in- 
teraction steps: the outer arm is produced by a first pas- 
sage of an intruder close to NGC 73 19 while the inne r 
plume formed in a more recent event (jMoles et al.lll997f) . 
This would explain that the second tail is still very young 
and not quite extended (regarding to the other one). 
Ku et all (|2005f ) pointed out that it took > 5 x 10 8 yr, 
after the close encounter with NGC 7319, for NGC 7320c 
to move to its current position. They argued that it is 
likely that the new tail is triggered by a close encounter 
between NGC 7319 and NGC 7318a, which is - 3 times 
closer to NGC 7319 than NGC 7320c is. 

The other important region is the pair of galaxies 
(NGC 7318a and NGC 7318b) on the west side of the 
group. The smallest one (NGC 7318b) has a veloc- 
ity of 5765 km s" 1 whereas the velocity of its counter- 
part (NGC 73 18a) is closer to th e r est of the group 
f6620 km s" 1 . ISulentic et~aTl l200l . IXu et al.1 (I2003L 
I2005D interpreted this difference as a shift of the first 
galaxy ~ 100 kpc toward us, after a head-on collision of 
the pair members ~ 10 7 — 10 8 yr ago. This scenario ex- 
plains the peculiar H I distribution outside the optical 
disk. Although the observational data do not permit to 
conclude on the origin of the filaments at the North of 
the pair (called NW tails, orange on Figure [1]), we note 
that these structures seem to expand in a cylindrical sym- 
metric way, as opposed to the usual tidal tails that are 
confined to a well-defined direction with respect to their 
host galaxy. This "ring-like" structure observed in SQ 
backups the hypothesis of a head-on collisi on between 
the two members of the pair (Higdon 1995). As men- 
tioned before, the recent collision of the intruder with 
the IGM triggered a large scale shock, between the pair 
and the main galaxy (purple on Figure [1} . 

Note also that the combination of internal and tidal 
processes shaped smaller pseudo-spiral arms in the main 
galaxy and the intruder (green and red on Figure [T]). 
The question of their origin is still to be address but it 
is likely that their formation has been triggered by tidal 
interactions. 

These tidal and gaseous structures observed today re- 
veal a complex gravitational history. Therefore, in ad- 
dition to observations, numerical simulations are a good 
way to test hypotheses and reproduce the evolution of 
such a complex system. The morphology and kinematics 
of the structures presented here and on Figure [1] are the 
main clues in our understanding of the quintet, and thus, 



the starting point of our modeling study. 

3. SERIAL APPROACH 

The most natural way to grasp such a complex ob- 
ject is to set global initial conditions. In this case, 
the simulation starts using an arbitrary time, since 
one can consider the progenitors as isolated galaxies. 
Thus, the game is to find the orbital parameters that 
lead to an interaction. This method has been widely 
used for the simulations of pairs or simple groups (e.g. 
i Toomre fc Toomrdll972tlBarneslll988l:lHibbard fc Mihosl 
1995HDubinski. Mihos. fc Hernquistlll996UBarnesll2004t 
iRenaud et al.l 12008). Indeed, it is rather easy to create 
individual models for each galaxy, to set the phase-space 
separations between them, and let them evolve, following 
the implemented laws of physics. However, dynamical 
friction and rapid evolution of the morphology add some 
difficulties to this exercise and a lot of attempts are often 
required before finding a good match with observational 
data. 

Another approach consists in exploring the parame- 
ters space with restr icted simulations and sele ct the best 
match among them ([Barnes &: Hibbard 12009). The very 
efficient simulation of point-mass galaxies surro unded by 
test particles (see e.g. lToomre &; Toomrelll972D allows to 
cover a broad range of parameters. iTheis h Kohld (|2001[ ) 
combined this method with genetic algorithms. Thanks 
to a maximum likelihood estimation, a good match with 
an observational data set can be rapidly found. How- 
ever, such a method faces problems to properly describe 
the orbital decay, an important characteristic of galaxy 
mergers, in particular for compact groups and thus, is 
not applicable here. 

In the classical A^-body approach, before interaction, 
it is still possible to approximate the trajectories of the 
centers of mass with a two-body solution. However, close 
passages would induce a loss of energy (dynamical fric- 
tion) and mass (creation of the tails) , that will invalidate 
the two-body result. Therefore, when a third galaxy en- 
ters the already affected system, its orbit (i.e. impact 
point and relative velocity) strongly differs from those 
one could expect analytically. That is why, for this study 
which implies complex trajectories, we use the traditional 
try-and-see method to reproduce the quintet. 

As SQ is currently at a typical intermediate stage of 
its evolution (and thus, displays a very complex web of 
features), it is very difficult to define a complete set of 
the numerical parameters which best describes the ob- 
served object. Therefore, we decided to follow a serial 
(as opposed to parallel) approach, assuming that most 
of the features of the quintet have been created sepa- 
rately and then, gathered without strong interaction or 
disturbance between them. This quite bold assumption 
allows us to split the problem into simpler interacting 
pairs, and to get a good idea of the historical evolution 
of compact groups in general. Our method first creates 
a galaxy-galaxy interaction and uses the result as one of 
the progenitor of the second step, and so-on. 

Luckily, SQ does lend itself to such a slicing in 
both spatial and t ime dimension s. Previous studies 
(jMoles et all Il997t IXu et all 12005ft suggested that the 
eastern tidal tails were built ~ 500 — 700 Myr ago for 
the longest one and ~ 200 Myr ago for the narrowest 
one, although the western pair of galaxies would have an 
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Fig. 2. — Schematic view of the major steps in the history of 
SQ. Very early phases and future (marked with ?) are subject to 
arbitrary assumptions that may not be realistic. See Section [5] 
however, for a discussion on the possible future of the group. 

interaction history of ~ 10 Myr only. 

Therefore, our serial approach consists in creating the 
two main tidal arms in the eastern region, involving 
NGC 7319 (the main galaxy), NGC 7320c (first collider) 
and likely NGC 7318a (second collider), and then add 
NGC 7318b (the intruder) in the western side for the last 
interaction event. Because of its separation with the rest 
of the group and the lack of linking tails, we assume that 
the first interaction of the remaining galaxy (NGC 7317) 
has not taken place yet. It may fall into the center of 
the group in the future but is not directly linked to it at 
present time. 

The interaction scenario, divided into six major steps, 
is schematically presented in Figured) (1) formation of 
the outer tidal tail, (2) creation of the inner (younger) 
tail, (3) western interaction between NGC 7318b and 
the rest of the group, (4) large scale shock, (5) involve- 
ment of the fifth galaxy (NGC 7317) and finally, (6) 
possible merger. According to this scheme, the cur- 
rently observed objects stand at the end of step (4). We 
do not have many clues about the last two steps but 
we can reasonably assume that SQ will fit in the gen- 
eral evolution t heory of compact groups, as described in 
Uohnson et al.l (|2007| ) , i.e. tidal harassment of its gaseous 
envelop and, possibly, formation of a massive ellipsoid 
from the progenitor s (see a discussion on th e shape of 
merger remnants in iWeil fe HeriicTuistl (l996l) . Depend- 
ing on the relative velocities of the galaxies, the final 
merger would involve all or only part of the members, and 
could occur in more th an a Hubble time, as suggested by 
Mamonl 119871 (see also [A thanas soula. Makino. fe Bosmal 
19971 and lAceves fc Velazquez] 120021 for a discussion on 
the overmerging issue of compact groups). It is likely 
that a similar approach applies to other compact groups, 
which undergo a sequence of collisions possibly leading 
to the hierarchical formation of a massive galaxy. 

From a simulation standpoint, the last step of gath- 
ering them into a single model is much more difficult to 
perform, even with a good description of the interactions 
of pairs of galaxies. Because the progenitors cannot be 
considered in an equilibrium stage after they have ex- 
perienced their first interaction, one has to carefully re- 
tune the initial parameters to get similar results than 
in isolation. NGC 7320c is already in a remote posi- 
tion when NGC 7318b enters the group, and thus should 
not be affected by the intruder. NGC 7318a however, 



stands very closely to the main galaxy. The trajectory of 
NGC 7318b must be re-set to compensate the attraction 
of NGC 7319 and keep the same impact point on the disk 
of NGC 7318a. In addition to this geometrical effect, one 
also has to consider the mass loss of the second collider 
from which a tidal arm has already been stripped. Dif- 
ferences with the simulation of the simple pair should be 
limited however, thanks to the high velocity of the in- 
truder: NGC 7318b undergoes the perturbations of dis- 
tant galaxies for a very short time, which reduces their 
effects. 

4. MODELLING STEPHAN'S QUINTET 
4.1. Numerical method 

For this study, we use a purely gravitational iV-body 
code to reproduce a global scheme of the interacting his- 
tory of SQ. Hydrodynamical approach will be added in a 
forthcoming worlfl. Each spiral galaxy is set up as self- 
gravitating disk + bul ge + dark matter halo model using 
the magalie software ([Boilv et al.ll2001[). included w ithin 
the Nemo stellar dynamics package (Teubcn 1995). We 
use extended bulge components, considering that they 
account for both the compact ce ntral bulge a nd the 
spheroid. The tree code gyrf alcON (Dehncn 2002) solves 
the equations of motion, with a Plummer softening ker- 
nel (e = 23 pc) and a tolerance parameter = 0.6. The 
conservation of energy has been verified by running a 
simulation of the model of NGC 7319 in isolation. We 
noted a relative difference of 10 -4 after 600 Myr. This 
small value gives us confidence in the results obtained 
with this code and this set of parameters. 

All our progenito rs, except NGC 73 20c, are made of an 
exponential disk, a iHernqui st ( 1990) bulge and a trun- 
cated isothermal sphere as dark matter halo. The first 
collider (NGC 7320c) does not show a clear internal tidal 
signature, hence we assume a spherically symmetric dis- 
tribution and only use the bulge and halo components. 
The mass ratio between all individual particles is set to 
unity, to avoid any two-body heating. The scaling of 
the galaxies in mass, velocity and length is in a such 
way that the average density and the virial parameter 
are conserved. Table [T] lists the parameters used for our 
best model. Considering that the initial mass (dark mat- 
ter included) of the main galaxy (NGC 7319) is about 
~ 2.5 x 10 12 M Q , our mass resolution is ~ 3 x 10 6 M . 
In the following, t = corresponds to the snapshot that 
best-matches the observations. 

4.2. Formation of the large scale tidal tails 

since iToomre fc Toomrd ((19721 ), it is well-known that 
the tidal features formed during an interaction are made 
of the material stripped from the galactic disks, because 
they are dynami cally cooler than the b ulges and the ha- 
los. H I maps by Will iams et all (|2002D suggest that the 
two major arms in SQ originate from the main galaxy 
(NGC 7319), as the greatest velocity difference between 
its disk and the tips of both arms is ~ 0.9% only. The 
elongation of the two arms clearly shows that the causes 
of their formation must be prograde (direct) interactions 
between collider(s) and a very large disk. 

4 We assume that the gas, which represents only a small fraction 
of the total mass, does not strongly affect the morphology of the 
system. 
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TABLE 1 

Parameters of the simulation. 



Parameter 


NGC 7319 


NGC 7320c 


NGC 7318a 


NGC 7318b 




(main galaxy) 


(1st collider) 


(2nd collider) 


(high-speed intruder) 


Number of particles 


870,000 


22,000 


300,000 


150,000 






Disk parameters 






Number of particles 


108,700 





37,500 


18,700 


Mass 


1.00 


0.00 


0.34 


0.18 


Vertical scalelength 


0.10 




0.07 


0.06 


Radial scalelength 


1.0 




0.7 


0.6 


Cut-off radius 


10.0 




7.0 


5.6 


Toomre parameter 


1.5 




1.5 


1.5 






Bulge parameters 






Number of particles 


54,400 


1,600 


18,700 


9,400 


Mass 


0.50 


0.10 


0.17 


0.09 


bcalelength 


2.0 


0.7 


1.4 


1.1 


Cut-off radius 


10.0 




7.0 


5.6 




D 


ark matter halo parameters 




Number of particles 


706,900 


20,400 


243,800 


121,900 


Mass 


6.50 


1.30 


2.21 


1.17 


Scalelength 


2.0 


0.7 


1.4 


1.1 


Cut-off radius 


20.0 




14.0 


11.2 






Disks inclinations 






(9x,9y) 


(0,0) 




(180,23) 


(0,-23) 






Initial coordinates 






(x,y,z) 


(0.0,0.0,0.0) 


(12.2,6.4,12.2) 


(-50.9,0.0,23.7) 


(-8.6,10.8,-133.0) 


(v x ,Vy,V Z ) 


(0.0,0.0,0.0) 


(-0.80,0.30,-0.50) 


(0.60,-0.06,-0.12) 


(0.04,0.00,1.00) 



Note. — The parameters are given in numerical units (G = 1). The conversion factors are 3.2 X 
10 11 M , 2.3 kpc, 775 km s" 1 and 2.9 Myr. The non-dimensiona l izatio n is based on the separation 
between NGC 7319 and NGC 7318a measured by Liscnfcld ct al. (2004) and the velocities given by 
ISulentic et"aT] |2001), assuming a recession velocity of 6400 km s 1 and Hq = 75 km s 1 Mpc 1 
which gives a distance of 85 Mpc. 



4.2.1. Outer tail 

The outer tail, which spans from the main galaxy to 
the south and then curves toward the east (see FigureH]), 
is obviously the oldest tidal structure in the quintet. In- 
deed, its very low surface brightness and large width sug- 
gest that it had time to expand and vanish into the IGM. 
Furthermore, the tidal strip ping of the associated H I 
distribution (noted Arc-S by iWilliams et al.l (|2002f) . see 
their Figure 5) by the collider implies that no major in- 
teraction had expelled the gas before, and thus that this 
event is the first in the history of the group. As the H I 
tail points toward the North, close to NGC 7320c, it is 
very likely that this galaxy is responsible for the first 
interaction (and thus justifies its nickname: the first col- 
lider). Both the H I distribution and the low brightness 
of the tail suggest that it has formed in a distant inter- 
action, during which the gas has been stripped off while 
only the less bound stars were tidally affected. 

In our modeling approach, this first step accounts for 
the creation of a long tidal tail originating from the main 
disk, a small counter-tail, and an orbit for the collider 
that drives it toward its current position. We note that, 
because of its mainly gaseous structure, the tail must 
count only few stellar particles of our model. To do so, 
we set NGC 7320c on a highly eccentric prograde orbit 
around NGC 7319, with a non-penetrating encounter, as 
shown in Figure [3J This setup has to bring it in the 
vicinity of the main disk but keep enough energy to send 
it away, at ~ 85 kpc from NGC 7319, as observed. 

There is no clear clue about the inclination of the or- 
bital plane with respect to the plane of the sky. The two 
progenitor galaxies show similar line of sight (LOS) veloc- 
ities, which suggests a weak (or even no) proper motion 




-40 -20 20 40 

x [kpc] 

Fig. 3.— Morphology of NGC 7319 and NGC 7320c after the 
first interaction (t = —260 Myr). The outer tail is very faint for 
our stellar model but should be much more visible with a gaseous 
component (similarly to the observations of H I). All the parame- 
ters used (including the orientation of the line of sight) are those of 
the final model of the entire quintet. The solid lines are the orbits 
of the center of mass of the two galaxies, with triangles marking 
their positions before the interaction. 

along this line, by contrast with NGC 7318b. Therefore, 
it is likely that both orbital planes are nearly identical. 
Note that the repartition of the measured LOS velocity 
between recession and proper motion of a given galaxy 
is unknown. If one assume a common group velocity 
(based on the one of the main galaxy), it seems rea- 
sonable to attribute the difference to both a geometrical 
offset along the LOS and a proper motion. In the case of 



6 



NGC 7319 and NGC 7320c, the velocity difference is too 
small (~ 30 km s _1 ) to allow any conclusion. Therefore, 
we first set the orbit of the first collider on the plane of 
the sky. 

In the Keplerian approximation (valid until the peri- 
center passage), the main galaxy stands at one of the 
two foci of the orbit of the collider. Because NGC 7319 
must undergo a prograde encounter, the direction of ro- 
tation of NGC 7320c along the orbit is determined by 
the spin of NGC 7319. Here again the choice is arbitrary 
because of the lack of observational clues. We set the ro- 
tation of the main galaxy to be counterclockwise, when 
viewed from the Earth, implying the same direction for 
the trajectory of the collider. Therefore, to bring it to its 
observed position, it must pass the main disk westward. 

4.2.2. Inner tail 

The material of the two tails clearly shares the same 
origin: the main disk. However, the inner tail shows 
a much higher surface brightness, is thiner and shorter, 
which suggests that it has been created after the outer 
one, and thus during a second interaction. A glance at 
the velocities and H I distribution of the members of SQ 
limits the possibilities to NGC 7320c and NGC 7318a. 

In the first case, NGC 7320c has to pass twice in the 
vicinity of the main galaxy: once to create the outer tail, 
and a second time to build in inner one. This constrains 
the possible trajectories to bound orbits, with enough 
energy to "escape" from the gravitational potential of 
main disk twice and finally go to the observed remote 
position (~ 85 kpc away from NGC 7319). However, in 
order to strip such a amount of material (H I and stel- 
lar tidal tails), the collider must pass very close to the 
center of the main disk, where the density is high. In 
these regions, the dynamical friction plays an important 
role and forbids the satellite galaxy to keep the same or- 
bital energy. Because the friction already makes it brake 
after the first pericenter passage, the orbital decay after 
two passages does not allow the required separation, as 
visible on Figure |4] To verify this, we have run a sys- 
tematic survey on the orbital parameters space, made 
of 96 simulations involving NGC 7319 and NGC 7320c. 
As no model gives the right final separation between the 
two progenitors after two passages, we conclude that it 
is very unlikely that the same galaxy creates both arms 
during two different interaction events. 

In another scenario, the collider and the main disk 
interact only once, and both tails are created at the 
same time. Indeed, the effect of the dynamical fric- 
tion is considerably reduced and the orbital decay be- 
comes less important than in the two passages scenario. 
With a highly inclined encou nter (an approach com- 
monly used for M 51, s ee e.g. iToomre fc Toomrd IT973I : 
Thcis & Spinncker 2003) the tidal bridge forms between 
the collider and the main disk, while a counter-tail is 
created in the opposite side of the disk, representing the 
inner and outer plumes, respectively. In our simulations, 
the mass ratio, the inclination of the orbital plane and 
the impact velocity have been adjusted to avoid the de- 
struction of the disk on the one hand, and the merger of 
the two progenitors on the other. After few tests, a good 
agreement with observations has been found for the po- 
sition and inclination of the tails, as shown in Figure [5l 
However this scenario fails to take into account the age 
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Fig. 4. — Hypothesis of the creation of two tails by a single col- 
lider. The figure plots the morphology of the system after the first 
pericenter passage. The outer tail is already well developed in the 
southern part. Solid and dotted lines shows the orbital decay in 
the trajectories of the collider and the main galaxy, respectively. 
Initial positions are marked with triangles. The second pericenter 
passage will lead to a strong loss of orbital energy and finally the 
merger of the two progenitors. Therefore, this scenario has been 
rejected. 
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Fig. 5. — Hypothesis of the creation of the two tails with one pas- 
sage only. Both structures have correct positions and inclination. 
However, the age difference between the plumes is not reproduced 
by this scenario. 

difference between the two plumes (because both are cre- 
ated at the same time) , and the a priori low inclination of 
the orbital plane with respect to the plane of the sky. In 
addition, both plumes present similar surface brightness 
when formed the same way, another difference with the 
observations, suggesting that this scenario is not correct. 

The last remaining possibility involves NGC 7318a. 
When considering this galaxy as the second collider, the 
age difference (~ 3 - 5 x 10 8 yr) between the two tails 
can easily be explained as a result of two well separated 
events. As the "two events - single collider" scenario re- 
veals itself to transfer too much of its orbital angular 
momentum through dynamical friction, we set two col- 
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Fig. 6. — Schematic view of the orbits of the two colliders around 
the main galaxy. According to the ages of the tails (and thus 
the ages of the tidal interactions) the two colliders might interact 
where their trajectories cross. Circles show the position of the three 
galaxies at present time. 



lider to hit the main disk sequentially. Because the outer 
tail is more extended and older, the cause of its creation 
had more time to escape the central region of the group. 
Therefore, NGC 7320c is likely responsible for it, while 
the inner tail would be created by NGC 7318a, few 10 s yr 
after. In this case, the outer tail is a bridge between the 
first collider and the main disk, while the inner tail is a 
classical arm (or counter-tail) from the interaction with 
the second collider. 

The numerical simulation of this scenario is more diffi- 
cult than usual, as the second interaction has to occur in 
an already disturbed system, without affecting the exist- 
ing structures (i.e. mainly the outer tail). Any Keplerian 
approximation (even before the pericenter passage) is not 
possible in this pseudo three-body problem. In addi- 
tion, the existence of extended tidal structures voids the 
point-mass approach. For simplicity however, we set the 
initial conditions with Keplerian orbits, keeping in mind 
that the actual trajectory will rapidly differ from those 
of our naive setup. In the previous section, we have set 
the spin of NGC 7319 to be counterclockwise. To get the 
prograde encounter for the main galaxy, the second col- 
lider must hit NGC 7319 southward before going toward 
the west. Unfortunately, the observational data reveal an 
extremely dense and complex web of gaseous and stellar 
structures between these two galaxies (recall Figure [T]). 
Hence, it is rather difficult to look for tidal signatures 
for NCG 7318a, that may indicate if its collision with 
the main disk is prograde or retrograde. However, a dif- 
fuse structure linking the two galaxies is visi ble in deep 
images (see e.g. the HST/WFC3 image by |NoH| l200l 
ER04) and thus, suggests the presence of a tidal bridge 
from NGC 7318a toward NGC 7319. Note that most of 
these tidal features have been strongly disturbed by the 
recent collision with NGC 7318b (see the next Section). 

According to their current positions, the two colliders 
have trajectories which could make them interact, be- 
tween the two collision events with the main galaxy. In- 
deed, after its collision, NGC 7320c flies eastward, north 
of the main disk, while NGC 7318a is approaching to cir- 
cle around the the main galaxy from north-east to south- 
west, as shown in Figure [5] The lack of observational 
clues (like tails between NGC 7318a and NGC 7320c or 
a common halo) suggests that the two galaxies have not 
interacted with each other. In order to keep the same 
timescales and the same orbits with respect to the main 
galaxy, the only solution to avoid the collision between 
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Fig. 7.— Morphology of NGC 7319 and NGC 7318a after the 
formation of the inner tail (t = —80 Myr). NGC 7320c is included 
in the simulation, but not shown here. The tail and the spread 
ISM between the two galaxies are clearly visible. A secondary tail 
originating from NGC 7318a is expanding toward the west but will 
vanish during the next steps of the interaction history of SQ. The 
outer tail, created by NGC 7320c, is already too difused to be 
visible. 



NGC 7318a and NGC 7320c is that they do not share 
the same orbital plane. According to the measured ra- 
dial velocities, the second collider might be behind the 
first one, at present time. This suggests to tilt the orbital 
plane of NGC 7319a so that it starts closer to the Earth, 
interacts with the main galaxy, and then goes behind. 
By doing this, the trajectories of the two colliders do not 
cross. 

This scenario involving three galaxies allows to tune 
the initial parameters in order to place the two intruders 
at their right positions after the interaction, i.e. in the 
far east for NGC 7320c and near the center, ready for 
a collision with NGC 7318b, in the case of NGC 7318a. 
The inclination of the orbital axis and the initial starting 
point along the trajectories have been chosen so that the 
two progenitors reach their observed position at the same 
time (i.e. now). The initial Keplerian orbit of the second 
collider has to be corrected to account for the mass loss 
and shift in position due to early collision. The final 
result is shown in Figure [7] 

We note the presence of a strong spiral and even a 
barred shape in NGC 7319 (see Figure [T]). Previous 
studies suggested th at such structures could resul t from 
a tidal interaction ([Walker. Mihos. fc Hernquistl 119961 : 
IBerentzen et al.1120041 ). On the one hand, it is very likely 
that more recent interactions erased or at least disturbed 
them. On the other hand, deep images clearly reveal that 
such asymmetric structures still exist in the quintet, in 
particular in the main disk. Our simulation managed 
to reproduce the spiral arms and even a weak bar af- 
ter the second collision. They are transient features that 
will vanish within a few dynamical times. By integrating 
NGC 7319 in isolation during 600 Myr, we have checked 
that these structures are not intrinsic to the model, but 
created by the interaction with the other galaxies. 

4.3. Formation of the western region 



For the next step in the interaction history of the quin- 
tet, NGC 7318b hits the triplet NGC 7318a-19-20c and 
its already existing web of tidal features. Although the 
velocity fiel d reveals itself to be extremely complex i n 
this region ()Sulentic et al.l 120011: IWilliaras etlfll2002j) . 
one can extract some hints from it. First of all, the 
large velocity difference measured between NGC 7318a 
and NGC 7318b (~ 850 km s" 1 ) clearly requires an or- 
bit highly inclined from the plane of the sky for the 
NGC 7318b (the intruder). Secondly, because the in- 
ternal tidal structures (spirals and bars) are still blurry, 
the interaction has taken place recently. Therefore, the 
velocity of NGC 7318b should be high enough to sepa- 
rate it from the rest of the group in such a small amount 
of time. Consequently, it is possible to consider the in- 
truder as a galaxy coming from the background at high 
velocity and hitting the group in a configuration close to 
the one currently observed. 

4.3.1. Tidal tails of the pair 

The most visible features in this region are two thin 
tails (NW tails in Figure [1} expanding northward from 
the two galaxies NGC 7318a/b. Their mean measured 
velocity of ~ 6000 km s _1 suggests that they stand be- 
tween the two galaxies. Therefore, these structures have 
been created by the tidal interaction of the second col- 
lider with the high speed intruder. It seems that each tail 
originates from one of two galaxies, the western one being 
from NGC 7318a while the other starts from NGC 7318b. 
Simulations of encounters show that a tidal structure ei- 
ther links the two galaxies (bridge) or point away from 
the counterpart (classical tail). Figure |S] shows the result 
of an interaction between two disks galaxies that leads to 
similar morphology than the pair. However, the northern 
tail of NGC 7318b is connected to its disk on the right- 
hand side, contrary to the observations of SQ. Hence the 
crossing of the two NW tails we see in SQ would be a pro- 
jection effect along the line of sight. In other words, the 
two tails would not share the same plane. However, this 
cro ssing matches t he position of the starburst region SQ- 
A (|Xu et al.lll999f) . It is likely that the over-density due 
to the physical overlap of the tails favors the formation of 
stars. Therefore, the crossing is not a simple projection 
effect but a volume common to both tails which have 
formed, in this case, in an unusual way. 

NGC 7318b shows another thin arm (SW in 
Figure [J) containing tidal dwarf galax ies can- 
didates (IHunsberger. Charlton, fc Zaritskvl 119961 : 
iMendes de Oliveira et al.l I2001D . The origin of this 
structure pointing toward NGC 7318a remains mysteri- 
ous, which confirms the complex history of this region 
of the group. 

A possible explanation consists in considering all these 
arms as a single ring-like structure. Face-on coll isions 
of disk galaxies create such rings (iHigdonl H99l and 
also explain the star formation along it ()Beirao et al.l 
I2009D. The velocity field measurements in this area 
(Will iams et~aT1 120021 ) place the ring between the two 
galaxies, but closer to NGC 7318b, suggesting that the 
external features come from its progenitor disk and un- 
derwent a significant brake regarding to the more com- 
pact nucleus. However, this scenario does not explain 
the asymmetry of this structure nor the crossing of two 
NW tails in the SQ-A area. For the moment, the obser- 
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Fig. 8. — Simulation of the western pair, in isolation. Here, the 
crossing of the two northern tails is real (i.e. not a projection 
effect). However, the main tail of NGC 7318b is not correctly 
positioned with respect to its nucleus. The crossing of the tails 
with a curvature toward the counterpart galaxy (and not away 
from it, as shown here), as observed in SQ, is very unlikely in a 
classical encounter. Hence these structures should have a different 
origin. 

vations reveal extremely complex morphology and veloc- 
ity field whose dynamical explanation seems still out-of- 
reach. However, it is certain that this region results from 
a recent high-velocity collision between a disk intruder 
coming from the background and NGC 7318a, which had 
already been altered by a previous interaction with the 
main galaxy. Consequently, one could expect the pro- 
genitor NGC 7318b to have a disk large enough to form 
the ring-like structure (or to contain all the material that 
has been spread away). 

Our goals in the modeling of this region remains mod- 
est, as we simply want to create two northern arms and 
possibly a third tail, southern. The fine details of the 
tidal structures, the crossing and the internal bar or spi- 
rals in NGC 7318b are out of the scope of the present 
paper. The results are visible in Figure [9l 

Tests showed that the required velocity corresponds 
to an unbound orbit of NGC 7318b with respect to the 
rest of the group, even when taking dynamical friction 
into account. Indeed, considering that the inclination 
of the disk of NGC 7318a is small with respect to the 
line of sight, the intruder undergoes a strong friction 
during only the short time it needs to go through this 
disk. In addition, the density of the IGM within which 
the intruder flies is low and hence, the dynamical fric- 
tion there does not play a crucial role. In this view, the 
high- velocity intruder does not loose enough of its orbital 
energy to remain in the gravitational well of the group, 
and escapes in the foreground. If this scenario is con- 
firmed, it might imply a different origin for the intruder 
than for the other galaxies in the group. Note however 
that these questions remain highly speculative and will 
not be addressed in this paper. 

4.3.2. Large scale shock 

The other major feature of the western part of SQ 
is the ~ 40 kpc long shock, between the high-speed in- 
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Fig. 9. — Morphology of the western pair shown in the plane of 
the sky (top panel) and after a rotation of 30° around the vertical 
axis (bottom panel). Again, the other galaxies (the first collider 
and the main disk) are included in the simulation but not shown 
here. The tidal structure on NGC 73f8b's disk exhibit a circular 
configuration, as opposed to a classical tail. This is well-visible on 
the bottom panel, with a spoke running southward. Such spoke 
could explain the SW-tail observed. 

trader and the main galaxy. Previous studies suggested 
that the shock results from the hig h-speed collision be- 
tween NGC 73 1 8b and the IGM dPietsch. et al.1 IT9971: 
Trinchieri et~aT1 120051: Ku et all 120031 : lAppleton et al.l 
20061) . In the scenario presented here, this medium comes 
from the interactions involving NGC 7319, NGC 7320c 
and NGC 7318a. It contains the gaseous and stellar 
components tidally stripped during the formation of the 
eastern tails and the "bridge" between the main galaxy 
and its second collider. Indeed, observations in sev- 
eral wavelengths reveal that this region is filled with 
atomic and molecular gas (see e.g. iLisenfeld et al.1 [20021 : 
lO'Sullivan et al.l[2009L and references therein). The re- 
cent intrusion of NGC 7318b in this medium triggered a 
hydrodynamic shock that propagates eastward. In addi- 
tion to the relatively high densities of the tidal tails, the 
shock might also help to trigger the star formation, as ob- 
served in SQ-A, via radiative transfer (see the suggestion 



of lBarneil2004 and his application to the Mice). 

Note that, because of its tidal origin, this gas stands 
between NGC 7319 and NGC 7318a (which is behind the 
main galaxy). Considering that the high-speed intruder 
comes from the background, it has hit NGC 7318a be- 
fore the gaseous medium (which is closer to us). In this 
respect, the shock is younger than the tidal features ob- 
served in the western region (NW and SW tails). It is 
clear that the origin, formation and propagation of the 
shock result from hydrodynamic processes, not discussed 
here. Hence, the shock is not modeled in the present 
work, but will be in forthcoming studies. 

4.4. Final model 

In Figure 1 101 we show the complete simulation of the 
group, at t = 0. Even if small-scale details are not 
perfectly reproduced, the large-scale configuration of the 
galaxies and the tidal features match well the observa- 
tions. We note that the bar and the spiral pattern in 
NGC 7319 are induced by external effects (the interac- 
tion with NGC 7320c and NGC 7318a) and not internal 
processes. Most of the IGM is stripped off the main 
disk but also originates from the late collision between 
NGC 7318a and NGC 7318b. 

Figure [11] displays the orbits of the three satellites in 
the reference frame of the main galaxy. Before its inter- 
action with the group, NGC 7318b has an almost radial 
motion and does not show a large displacement in the 
plane of sky. However, after the collision, the attraction 
of the two massive galaxies (NGC 7319 and NGC 7318a) 
deviates the intruder from its initial orbit. This explains 
the peculiar shape of its trajectory, in projection. 

As in the observations, the outer plume is very faint 
(and almost not visible) in our simulations. Still, the 
external part of the main disk has been spread toward 
NGC 7320c, in a very diffuse fashion. Almost all the stars 
of this tail visible in Figure [TU] are between NGC 7320c 
and the tip of the inner plume. A comparison with H I 
observations suggests indeed that this tail would be made 
of mainly the gaseous envelop of the main galaxy and 
thus, would not count a lot of stars. New runs includ- 
ing hydro dynamics are curr ently performed to confirm 
this idea (Hwang et al. 2010). Note that for more of the 
stellar component to be ejected, the impact parameter 
should be smaller, but the main galaxy would be affected 
too strongly to permit the creation of the second, inner 
tail. 

Our model fails to place the western pair of galaxies 
at its observed position, i.e. southern than in our sim- 
ulation. Note that this problem cannot be corrected by 
simply changing the orientation of the orbits, because 
this will imply a rotation of the tidal tails too. An incli- 
nation of the orbital plane, leading to projection effects 
could be a solution, although no observational element 
indicates such a configuration. 

The fine structures in the west pair are not well re- 
produced. Indeed, a description of such features requires 
a very precise setting of the parameters. Unfortunately, 
NGC 7318a is not in equilibrium (after its encounter with 
the main galaxy) when it hits the intruder, as many tran- 
sient features have not dissolved yet. Therefore the im- 
pact point can not be precisely defined, nor the relative 
velocity or the mass ratio. In addition, most of the struc- 
tures observed seems to have an important (even domi- 
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Fig. 10. — Morphology of the simulated Stephan's Quintet, at t = 0. This figure stacks the previous Figures \3\ [7] and [9] this time showing 
all the galaxies. The additional galaxies NGC 7317 (west) and NGC 7320 (foreground) have been manually overlaid to the configuration 
of the rest of group, for graphical purposes only. Note that these two galaxies are not included in the simulation itself, because of lack of 
observational signatures of past interactions that may constrain their trajectories and velocities. 
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Fig. 11. — Orbits of the first collider (dashed), second collider 
(solid) and the high-speed intruder (dotted) in the plane of the 
sky, with respect to the position of the main galaxy. Initial posi- 
tions are marked with triangles. The arrow indicates the direction 
of NGC 7318a, whose initial position is out of the frame of the 
Figure. The change of curvature on the orbit of NGC 7320c re- 
flect the modification of the central potential of the group, due to 
interactions. The merger of NGC 7318a is clearly visible after its 
second pericenter passage. 



nant) hydrodynamic part, that cannot be reproduced in 
purely gravitational simulations. The asymmetry of the 
system and the over-densities matching the positions of 
the NW-tails are correctly reproduce however, thanks to 
the formation of the ring-like structure. 

In addition to the morphology, one has to check that 
the kinematics of the group is also well-reproduced. Fig- 
ure displays the velocity along the line of sight for 
the galaxies included in the simulation. If one assumes 
that the kinematics of the gas follows that of the stars, 
it becomes possible to compare the velocity patterns ex- 
isting in our model to those measured from H I data in 
previous studies. All the features show velocities close 
to the ones measured, in particular the gradient between 
the high-speed intruder and the second collider. The ve- 
locity p attern of the NW and SW tails also matches the 
data of iWilliams et al.l (|2002l . their Figures 10, 11 and 
12). 

5. DISCUSSION 

5.1. Assumptions of this model 

Despite a number of observational data about SQ, we 
have made some assumptions for the simulation pre- 
sented above. First of all, the distribution of dark matter 
is not well constrained by the observations. Here, we con- 
sider one spherical halo per galaxy. Note however, that a 
more complex di stribution involving triaxial cuspy g alac- 
tic halos (see e.g lHavashi. Navarro, fc Spr ingcl 20013) an d 
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Fig. 12. — Velocity along the line of sight, of the final model. A systemic velocity of 6650 km s _1 (yellow) has been assumed. The shift 
of the high-speed intruder with respect to the rest of the group, and in particular NGC 7318a is clearly visible. The ring-like strucutre also 
contrasts well with it enviroment. 



even a global group halo is possible. Changing the dark 
matter distribution in the simulation will require a fine 
tuning of the orbital parameters of the galaxies because 
the halos, being the most massive and gravitationally 
dominant components of the progenitors, set the details 
of their dynamical evolution. However, this would not 
affect the chain of events described above, nor the role 
of each galaxy. 

Similarly, we have assumed that the rotation of 
NGC 7319 is counterclockwise. Because its disk lies in 
a plane close to those of the sky, the rotation curve is 
not measurable, which makes our hypothesis purely ar- 
bitrary. Indeed, if one considers a clockwise rotation in- 
stead, the orbits of the two colliders simply need to be 
flipped, in order to keep prograde encounters. In this 
case, NGC 7320c passes southward of the main galaxy, 
while NGC 7318a goes northward (corresponding to a 
180° flip of Figure [6] around its horizontal axis). By sim- 
ply changing the initial conditions, it is possible to reach 
a comparable final configuration, which shows that the 
model presented here is not unique. The sequence of 
events and the origins of the structures are, however, 
very likely to follow our dynamical scenario. 

Finally, the mass ratio between the progenitors is not 
strongly constrained. In our model, we set it thanks to 
the relative luminosities of the galaxies (assuming a con- 
stant mass-to-light ratio) and then tuned it according to 
the dynamical results obtained during our exploration of 
the parameters space. One can reasonably assume that 



NGC 7318a and NGC 7318b, as observed today, have 
similar masses (ratio close to 1:1). Taking into account a 
mass loss of 30% of the second collider during its inter- 
action with the main disk, one obtains that NGC 7318a 
must be initially more massive than the high-speed in- 
truder. 

5.2. Future of SQ 

The sequential splitting of events matches the spatial 
distribution of the galaxies, in the sense that interactions 
are well separated in time and space. Figure [T3] plots the 
evolution of the distance between each galaxy and those 
of NGC 7319. This clearly shows that the interaction 
events (i.e. the low distances) follow a sequence in time. 
It also allows to estimate the age of structures like, e.g. 
the large-scale shock, once a triggering event is assumed, 
by searching for, e.g. a minimum distance between the 
progenitors. We note that our values (~ 150 Myr for the 
age of the outer tail, ~ 80 Myr for the young one, and 
about 10 Myr for the shock) match well the estimations 
from the literature presented earlier. 

Figure [13] shows that NGC 7318a merges with the main 
galaxy at t « 100 Myr. Simple estimates of escape ve- 
locities reveal that the mass of all the simulated galaxies 
but the intruder (3 1012 M) corresponds to an escape 
velocity of ~ 250 km s _1 at a distance of 200 kpc, i.e. 
less than the estimated velocity of the intruder. This 
is confirmed by our simulation: NGC 7318b is unbound 
with respect to the rest of the group. Note that Fig- 



12 



ro 100 



a 
z 











- - — : 


: \ / 

1 A 
' ' \ 












shock > 


A 


NGC 
NGC 
NGC 


7320c : 

7318a ; 

7318b 







200 
time [Myr] 




400 



Fig. 13. — Evolution of the 3D distance between SQ members and 
its main galaxy NGC 7319, in log-scale. Dashed, solid and dotted 
lines correspond to NGC 7320c, NGC 7318a and NGC 7318b re- 
spectively. This figure shows the sequential scenario of the creation 
of SQ, as all the episodes of the interaction history are clearly visi- 
ble and well-separated in time. The arrow marks the approximative 
time of the shock between NGC 7318b and the IGM. Again, t = 
corresponds to the best-match to the observations. 

ure [13] suggests that this is the case of NGC 7320c too. 
However, the remaining galaxy NGC 7317 (not simulated 
here) does not yield a high LOS velocity and is close to 
the central region of the group. Therefore, it will likely 
merge with the main galaxy and thus modify the gravita- 
tional potential of the group. This may change the fate of 
the escaping galaxies, by slowing them down and maybe 
making them fall back into the group core. Following our 
estimates, a central mass of ~ 2 x 10 13 M Q is required 
to bind NGC 7318b which yields a LOS relative velocity 
of 850 km s™ 1 in both the observations and our model. 
Unfortunately, the lack of tidal signatures and estimates 
for the tangential velocity of the missing galaxy forbids 
us to go further than speculations. 

6. CONCLUSION AND PERSPECTIVES 

Using collisionless TV-body simulations, we studied 
some possible formation scenarios of the Stephan's Quin- 
tet. We have run more than 3,000 attempts to simulate a 
fraction or the entire quintet and constructed a dynam- 
ical scenario based on the separation in time and space 
of the different interaction events. This scenario can be 
summarized as follow: 

• mild interaction between the first collider 
(NGC 7320c) and the main disk (NGC 7319). The 
old, faint tidal tail is created. Because it is difficult 
to see in our simulation, it may mainly consist of 
the H I external envelop of the main galaxy. 

• second interaction, between NGC 7318a and 
NGC 7319. The younger tail is built and points 
away from the second collider. A significant mass 
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fraction (~ 30% of the disks) of both galaxies is ex- 
pelled into the IGM. Internal structures are created 
in the main disk. 

• frontal collision of the high-speed intruder 
(NGC 7318b), coming from the background, with 
the second collider (NGC 7318a). Several arm- like 
tidal features are created, possibly via the forma- 
tion of a ring-like structure. 

• interaction between the high-speed intruder and 
the IGM. A large scale shock front forms. 

• interaction between NGC 7317 and the rest of he 
group. This event has not taken place yet and is 
thus not implemented in our simulations. 

After an exhaustive exploration of the parameters 
space, it appeared that other hypotheses for the forma- 
tion of the tidal tails are very unlikely. In particular, the 
origin and relative age of the outer and inner tails seem 
to be very well constrained. The case of the western 
structures should be investigated further when including 
hydrodynamics aspects. 

Such a scenario clearly shows a splitting of interac- 
tion events in space and time, over ~ 100 kpc and a few 
~ 10 8 yr, a property which may characterize the evolu- 
tion of compact groups. In the special case of SQ, the 
configuration observed reveals past episodes, whose im- 
prints are the tidal tails, current events like the collision 
of the western pair of galaxies, and future interactions 
e.g. between NGC 7317 and the rest of the group. In 
slightly more evolved groups, each event would be char- 
acterized by a tidal signature whose age may denote this 
time sequence. However, at the end of their hierarchical 
evolution, the light of compact group is mostly domi- 
nated by the central galaxy, which makes diffuse plumes 
more difficult to detect. 

Future numerical studies will explore the response of 
a g aseous component to the multiple interactions of 
SQ (jHwang et alJ l20Toh . This might bring new clues 
about the formation and early evolution of the large- 
scale shock. With a detailed hydrodynamical model cou- 
pled to our gravitational description, we could be able to 
study the large-scale shock, the starburst regions, and 
furthermore, to answer a lot of questions raised by a 
growing number of observations over a broad range of 
wavelengths. 
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